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We show that at sufficiently large chemical potential SU{N) lattice gauge theories in the 
strong coupling limit with staggered fermions are in a chirally symmetric phase. The proof 
employs a polymer cluster expansion which exploits the anisotropy between timelike and 
spacelike directions in the presence of a quark chemical potential fi. The expansion is shown 
to converge in the infinite volume limit at any temperature for sufficiently large /i. All expec- 
Qh' tations of chirally non-invariant local fermion operators vanish identically, or, equivalently, 

their correlations cluster exponentially, within the expansion. The expansion itself may serve 
as a computational tool at large n and strong coupling. 



^e-mail: tomboulis@physics.ucla.edu 



1 Introduction 

Elucidating the rich structure of the QCD phase diagram, and more generally that of SU{Nc) 
gauge theories, as a function of temperature and finite density continues to present a con- 
ceptual and technical challenge. A major obstacle in this endeavor has been the inability to 
perform lattice gauge theory ab initio simulations by standard Monte Carlo methods since the 
fermion determinant becomes complex in the presence of finite baryon chemical potential, the 
infamous "sign problem" . Significant progress has been made in recent years in overcoming 
some of these problems by a combination of numerical simulations, mostly for small chemical 
potential, analytical techniques and investigation of model systems (see [1] [2] for review). 
Despite such progress, however, away from the region along the temperature axis the QCD 
phase diagram remains, for the most part, conjectural. In the absence of results established 
from first principles by simulations it is important to attempt to obtain any first principle 
results by other means. 

In this paper we consider SU{Nc) lattice gauge theory with Nf fiavors of massless stag- 
gered fermions in the strong coupling limit at temperature T and quark chemical potential 
fi. The theory has been investigated in the literature, in particular for the cases of SU{2) 
and 5C/(3) with Nf = 1, using a variety of approaches. Integrating out the gauge field at 
strong coupling leads to a representation of the partition function in terms of monomers, 
dimers and baryon loops [3], or monomers, dimers and polymers [4], [5]. In this representa- 
tion the sign problem is partly evaded allowing simulations. In such simulations in the case 
oi Nc = 3 ^ a clear sign of a chiral symmetry restoring first order transition was found at 
some critical ji. Similarly, for Nc = 2, restoration of chiral symmetry at large fx and/or T was 
seen in [5]. More recently, the two-color {Nc = 2), Nf = 1 case was investigated in [6] using 
the dimer-baryon loop representation with a new updating algorithm [7]. A second order 
transition to a chirally symmetric phase at some critical fi was seen in good agreement with 
mean field predictions. In the case of A'c = 3 such improved simulations were carried out in 
[8]. Another line of investigations follows a mean field approach. In [9], |10) . |11| an effective 
action was obtained by performing a 1/d expansion in the spatial directions and retaining 
only the leading terms, while leaving the timelike directions intact. Introducing auxiliary 
bosonization sources for condensates of interest, and finally integrating out the fermion and 
timelike gauge fields, treating condensates in a mean field approximation, leads to an effective 
action which allows a detailed study of the phase diagram in T, /i and quark mass m. This 
phase diagram exhibits a chiral phase at large /i and/or T. 

Here we will obtain an actual general proof of the existence of a chirally symmetric phase 
for the SU{Nc) gauge theory with Nf fiavors at strong coupling at chemical potential above 
a critical value. The basic approach adopted is as follows. Nonvanishing chemical potential 
and/or temperature introduce in the action an anisotropy between timelike and spacelike 
directions. This anisotropy grows with increasing //, and T, and provides the basis for setting 
up a systematic expansion in which the spacelike part of the action is expanded in a fermion 
hopping expansion in the measure provide by the timelike part. 

This type of expansion was first used in [12] to rigorously prove the restoration of chiral 
symmetry in lattice gauge theories at large T (and /u = 0) for all values of the gauge coupling. 
Here we exploit the anisotropy to set up an expansion in the presence of nonvanishing /x. For 
technical reasons [13], however, the expansion has to be set up in a somewhat different way 



than in [12]. We only consider the strong coupling limit. (This restriction is actually related 
to the complex nature of fermion determinants.) The resulting expansion can be formulated 
as a polymer-type cluster expansion |14j . |15j . The standard convergence criteria for polymer 
expansions can then be applied. One thus finds that the expansion converges in the infinite 
volume limit in spatial dimension d > 1 for chemical potential /i large enough. Within its 
region of convergence the expansion describes a chirally symmetric phase. 

In this paper we are only concerned with establishing the convergence of the expansion 
and thus the existence of a chirally symmetric phase at large chemical potential. In fact, 
the actual bounds we obtain (section 2.3 below) on the convergence region are rather non- 
optimal; they can be easily improved. We will not consider here the systematic evaluation 
of graphs in the expansion and/or resummation of classes of graphs, i.e. application of the 
expansion as a computational tool. It is interesting to note in this connection that quantitative 
results obtained by retaining a leading approximation in our expansion can be expected to 
be similar to those in [9], [10], [11] where only the leading terms in a 1/d expansion in the 
spatial directions [16] are retained. 

The paper is organized as follows. In section 2, after some preliminaries, the expansion 
is set up and its diagram structure laid out. We then formulate it as a polymer expansion 
of the logarithm of the partition function and of expectations of observables (section 2.2). 
Its convergence is examined in section 2.3. The preservation of chiral symmetry within 
the expansion follows then as a standard consequence of convergence of a cluster expansion 
(section 3). Some discussion and outlook for further work are given in section 4. The 
derivation of some explicit formulas used in the text is provided in the Appendix. 

2 The cluster expansion 

2.1 Preliminaries 

The lattice A is a ((i+l)-dimensional periodic hypercubic lattice of size Lf x L. It is convenient 
to always take the lattice lengths in Euclidean time (L) and space (Lg) to be even. A may 
be taken anisotropic with timelike and spacelike lattice spacings a,- and a^, respectively. 
This is the standard way of allowing independent variations of couplings and the physical 
temperature which is defined as T = (La,-)"^. The spatial lattice obtained as a particular 
equal time slice of A will be denoted by A^ . 

We use standard lattice gauge theory notations and conventions. Lattice site coordinates 
will be denoted by x = (x-^) = {x^, x) with A = 0, 1, . . . , d, and x = (x-^), j = 1, . . . ,d. We 
also often write x^ = t. Lattice unit vectors in the space and time directions will be denoted 
by j and 0, respectively. We generically denote lattice bonds by b, plaquettes by p, etc. Bonds 
will be specified more explicitly as 6 = (x, A), or bg = {x,j) =< x,x + j > if spacelike, and 
br = (x,0) =< x,x + > if timelike. Correspondingly, the gauge field variables Ub defined 
on each 6 G A will, as usual, often be more explicitly specified by Uj{x) and Uo{x). 

The gauge fields Ub € Gc are elements of the gauge group (color) Gc and taken to trans- 
form in the fundamental representation. Here Gc = SU{Nc). Fermions are introduced by 
associating on each site 21/ generators '4>a{x), tpa{x), (a = 1, . . . , z^), of a Grassmann algebra. 
The fermions are taken to transform as Nf copies (flavors) of the fundamental representation 



of the gauge group. In explicit enumerations we use a = l,...,A'^c)^ = l)---) ^f for color and 
flavors indices. We assume periodic boundary conditions for the gauge field and antiperiodic 
boundary conditions for the fermions. 

We consider fermions in the presence of chemical potential fi/ur- The lattice action for 
massless fermions is then given by: 

Sf = ^ o{ — ) ^ixhj{x)Uj{x)ip{x+j)-ijj{x + j)-fj{x)U^j{x)7p{x) 

+ Yl l[i'(.x)lo{x)e^'Uo{x)'^lJ{x + ^)-iJ{x + 0)Jo{x)e-^'U^Q{x)iJ{x)Y (2.1) 

br = {x,0) 

The matrices 7(6] = 'yx{x) defined on each bond b = (x, A) satisfy HfeeapTi^] — ^ fo'^ each 
plaquette p. For staggered fermions 

j^{x) = {-l)^^<^-\ 7o(x) = l. (2.2) 

For staggered fermions then one has z^ = NcNj, and thus a = (a, i). Recall that Nf staggered 
flavors correspond to ANj continuum flavors. For naive fermions 'yx{x) = 7^ are elements of 
the Euclidean Dirac-Clifford algebra satisfying {7^,7'^} = 26^'^! and are hermitian matrices. 
Naive fermions are of course unitarily equivalent to 2^^^'''^ copies of staggered fermions. We 
only consider staggered fermions in the following even though we often state formulas for 
general 7[6]'s. 

With Nf flavors of staggered fermions the fermion action in ()2.ip is invariant under a 
U{Nf) X U{Nf) global symmetry corresponding to independent rotations of fermions on even 
and odd sublattices [T7]. Explicitly, for any element {u,v) E U{Nf) x U{Nf), they are given 
by iIj{x) -^ uip{x) and 'ip{x) -^ 'ip{x)v'^ for even sites and ip{x) -^ vil>{x) and '4>{x) — > '4){x)u^ 
for odd sites. This symmetry is referred to as chiral symmetry. 

We work in the strong coupling limit, i..e. at inverse gauge coupling /? = 0, so a gauge 
field plaquette action term is absent. The full measure is then given by 

d/iA = Zl^ JJ dUb W di){x)dil){x) expi-Sp) , (2.3) 

beA xeA 

where, as usual, dUb denotes normalized Haar measure on the group Gc, and d'0(x)(iV'(x) = 
Y\^d'il}a{x)d'ipa(x) is the standard measure on a Grassmann algebra. The partition function 
Z/^ is defined by f dfj,\ = 1. The measure (12. 3p is invariant under the global chiral symmetry 
(as well as, of course, under the local Gc gauge symmetry). 

Expectations of general fermionic observables O are then given by 



0) = j dfiAO. (2.4) 

We are interested, in particular, in the possible formation of condensates of products of local 
fermion operators such as the chiral order parameter ip{x)ip{x) and various diquark or multi- 

quark operators. Introducing the notation x{x) = ( j ( \T 1 ' such operators are of the 

general form 

o = n x{xyt.,x{x) = n ^(^) ' (2-5) 

x&O x£.0 



where O (the support of O) is a finite set of sites, and T^ is some (unitary) matrix. 

2.2 The expansion 

The presence of a nonvanishing chemical potential in ()2.ip introduces an anisotropy between 
the spacelike and timelike directions which can be exploited, in a manner analogous to that 
of the case of large T [12], to set up a convergent expansion for large /i. For this purpose we 
rewrite (|2.ip more concisely in the form 

Sf = ^V5(x)Mg(C/)V'(2/) +5^Vi(x)Mg([/)]V(2/) , (2.6) 

x,y x,y 

where the matrices M'*-* and M*^*' have nonvanishing elements only between nearest-neighbor 
sites, i.e., on bonds, given by 

^%+3)iU) = \{ar/ash,{x)U,[x) and MJJ^.^JC/) = -i(a,/a,)7,(x)C/j(x) (2.7) 
for spacelike neighbors, and 

Mx,U(^) = ^7oe^f/o(x) and M^^^JC/) = -]^^,e-^ul{x) (2.8) 

for timelike neighbors. 

To set up our expansion we rewrite the measure ()2.3|) as follows. For each spacelike bond 



hg =< x,x -\- j > and each a let 

/i" = M^){Mll^^^i;{x + j)), , /J° ^ {i;{x + j)m[1\^^J^Mx) ■ (2.9) 

Note that {fl'_^ ")2 = 0. One then has 

2 u 
1=1 a=l 



exp(^(x)M(^()^^.)(^)^(x + i) + V;(:E + j)M[J^.)^(C/)V'(x)) = 1111(1 + /^") • (2-10) 



For each site x in a fixed time slice A^ define the measure 

(i/xx = ^n(iC/o(T,x)#(r,x)dV'(r,x)exp K^ V5(r,x)Mg^)(^, ,^)(C/o)V(t',x) (2.11) 

^ -r=l \t,t' ' ' J 



Note that 



J]V5(r,x)M[^;^)(^,^^)([/o)V(r',x) 



T,T' 

L 



J^[V5(r,x)7oe^[/o(r,x)^(r + l,x) -^(r + l,x)7oe-'^C/o^(T,x)V(T,x) (2.12) 



T=l 



with -(/^(L + ljx) = —ip{l,x), ^(L+l,x) = — ^(l,x) due to the antiperiodic ferniion boundary 
conditions. The factor z, defined by Jd/^x = 1, is then a partition function along a 1- 
dimensional timelike fermion chain given by 

L 

^= / n^^o(T,x)DetMi,*)([/o) (2.13) 

T = l 

with Mx (C/q) = \'^(t:x.)(t' x)(^o) I denoting the restriction of NVxy{Uo) to the submatrix of 
timelike bonds at fixed x. 

Using ()2.10p and (j2.1ip the full measure ()2.3p can now be expressed in the form 

2 1/ 

d^,^ = z-^^ n ^-"^ n '^^^^ n n n (1 + ^'o (^.u) 

xGAa fesSA 6s6Ai=lo=l 

with 

Za = Za/zI^^I . (2.15) 

The timelike part of the measure (|2.14p factorizes in a product with each factor representing 
the fermionic degrees of freedom coupled in a 1-dimensional timelike chain at fixed spatial 
coordinates x. It is very convenient then to adopt the gauge, frequently referred to as 
Polyakov gauge, where the bond variables Uq{t,-x) are independent of r and diagonal: 

C/o(r,x) = diag(e*^i(^)/-^,e^^2(x)/L^... ^gi0jv,(x)/L) _ (2.I6) 

= exp(ze(x)/L) (2.17) 

with 

AT, 

Y,Oa{^)=^ and e(x)=diag(0i(x),02(x),--- ,e7V.(x)). (2.18) 

a=l 

This is always possible by gauge transformations setting the C/q's at fixed spatial coordinates 
X equal and rotating into the Cartan subgroup [18j. After a time Fourier transform the 
timelike action (j2.12p becomes 

^V'(7",x)m|^|^^^^,^^^(C/o)V'(t',x) = ^ ^(/cm,x)i7osin(^A;„ + e(x)/L-i/_fjV'(/i:m,x), 

ry fcmSBZ 

(2.19) 
and, hence, the timelike propagator (covariance) in the background of the gauge field (I2.17P 
is given by 

a,,K0(x))^C(r-T',G(x)) = i Yl e'^"^(^-^'h{k„^,Q{^)) (2.20) 

fc^eBZ 

with 

C{k,e{ji.))=\i-fosm(k + e{jc)/L-in)Y . (2.21) 



In (|2.19p and (12.20p the summation is over momenta in the Brillouin zone (BZ): 

km = (2m - l)7r/L , -L/2 + 1 < m < L/2 ^ -tt + vr/L < /c^ < tt - vr/L (2.22) 



(integer m). Evaluation of (|2.20|) (Appendix A) gives 

/ /N p-iOa{-yi){T -t')/L . ,. 



and 



for (r - r') > 0, ^ > (2.23) 



C(r-T',e(x))„,,, = -<5,,<5.,[l-(-l)l^-^l] ,. .^ . ^-l.[L-\r-r\] 

for (r - r') < 0, ^ > . (2.24) 

For ^ < (i.e., nonvanishing antiquark chemical potential) replace fi,6a{x) by |/u|,— 0a(x), 
and reverse the sign condition on (r — r') in ()2.23p - ()2.24p . 

Note that C{t, 0(x)) vanishes for even r. This is a consequence of the chiral invariance of 
the action. The other salient property of C{t, S(x)) is its exponential decay for nonvanishing 
/i. This decay is in fact the crucial property needed for the convergence of the expansion 
below. It is worth pointing out that it can be extracted rather simply from ()2.20p . without 
its actual evaluation, as follows. Consider the quantity 

C{t - t', e(x)) = r — e^^-(^ - ^')c{k, G(x)) (2.25) 

J~7T 2vr 

obtained by replacing the Brillouin zone in p.20p with the interval [— 7r,7r], i.e., its infinite 
L (zero temperature) limit. Now, C{k,Q{x.)) is analytic in k for |ImA;| < |^| and bounded 
uniformly in Re A: G [— 7r,7r]. The periodicity of the integrand then allows one to "lift" the 
contour of integration in (j2.25p (a standard trick) from the real axis to, say, Imk = ^IfJ-l or 
— jIa*!, which immediately establishes the exponential fall-off of (|2.25p and the bound 



|(7(T,G(x))| < (sinh(|^|/2))-iexp(-^|r|) (2.26) 



It 
2 

in all cases. The exponential decay of C{t, 0(x)) then follows from the relation (see (|A.lip ): 

oo 

C{t, e(x)) = ^(-l)«C'(r + nL, e(x)) . (2.27) 

n=0 

Explicit evaluation, i.e., eqs. (|2.23p . (j2.24p . results in a better bound, but, if all one is 
interested in is demonstrating convergence at sufficiently large /i, the bound (j2.26p is sufficient. 
We also note (cf. Appendix A) that evaluation of (j2.25p and use of (j2.27p allows an alternative, 
simpler derivation of (|2.23p - (|2.24p than the direct evaluation of (|2.20p . 

We also need the quantity DetMx (?7o) resulting from integration of the fermions along 
the timelike chain at fixed spatial coordinates x in (|2.12p . Its straightforward evaluation in 
the gauge (|2.17p (Appendix A) gives: 

DetMi*^([/o) = DetC7~i(G(x)) 






2-'^^e 




(2.28) 



Note that in the case of A'c = 2 one has 9i = —02 = 9 and (j2.28p becomes 

DetC~\e{x)) = 2""^^-^^ [cos 9 + cosh Lfif^f , (2.29) 

which is indeed real and non-negative. 

Our expansion is generated by expanding the products over spacehke bonds in (|2.14p 
around the measure provided by the timehke part in (|2.14p . In other words, we will perform 
a hopping expansion in the spacelike directions with the spacelike hops connected in the 
timelike directions by the propagators (|2.23p - ()2.24p . In the following we set ar/cLg = 1 for 
convenience. 

Expanding first the products over the spacelike bonds gives 

B •' xGAs hs£B bs,l,a£B 

The sum is over all subsets B of the set of 2z^ copies of the set of space-like bonds in A. Each 
such B may be decomposed into a number of connected components where connectivity is 
defined as follows: two elements of a set B are connected if they may be joined by a sequence 
of time-like bonds. Corresponding to this decomposition of each set B, and after carrying 
out fermion and gauge field integrations, each term in (|2.3Up decomposes into a product or a 
sum of products J19j of factors, each factor being the value of a connected diagram. 

A connected diagram 7 consists of a number of space-like bonds connected by 1-dimensional 
chains of time-like bonds. It is important that the number of /^'" residing on each spacelike 
bond bs in a given diagram must result in an equal mod Nc number of Uj,^ and U^ 's on bs 
in order to obtain a nonvanishing result upon performing the Ub^ integrations. This is, of 
course, a characteristic feature of any fermion hopping expansion at strong coupling (absence 
of plaquette term) . Examples are shown in Fig. [TJ 

Each such connected diagram 7 will be termed a polymer. The value of the diagram Cil) 
is the activity of the polymer 7: 



C(7) = fU^^'-Ild^^s n fbT (2-31) 

= (H^^')/ n^^o(0W) nDetc"'(0W)n n a..,(e(x))/,.(2.32) 

xe7 xe7 xg7 ^^1 (TiTj) 

Here {x G 7} denotes the set of spatial sites obtained by projecting 7 onto a spacelike slice 
As (cf. Fig. [2]). In (j2.32p (TiTj) denotes pairs of points along a timelike bond chain at x 
connected by propagators, and I^ is a Nc, A^/-dependent constant resulting from the fermion 
and Ubg integrations in (I2.3ip . The expansion (j2.30p then becomes: 

00 ^ k 

^a = 1 + Ef ^ nC(70, (2.33) 

k=l {7i,-,7fc)6X'fe i=l 

where V^ denotes the set of all sets (71,72, ••• ,7fc) of k disjoint polymers. In (j2.33p we 
conveniently sum over all ordered sequences of polymers and, correspondingly, divide by kl. 
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(a) 



(b) 



(c) 



Figure 1: Examples (a) - (c) of connected diagrams built of spacelike bonds (solid lines) 
connected by timelike propagators (broken lines). An equal mod N^. number of U and C/^ 
(indicated by arrows) must occur on every spacelike bond [N^ = 3 in (b)). Note that diagrams 
such as the diagram in (c) cannot be factored in the product of two diagrams (a) since 
propagators along a common set of timelike bonds are "welded" together by the integration 
over the background gauge field residing on them (cf. ()2.32p ). 



The log of Za is now obtained through an application of the moment-cumulant formalism 
of statistical mechanics. Let X = (71, ... , 7,^.) be a set of k polymers (not necessary distinct); 
we will refer to such as set as a /c-polymer (fc-cluster). To each such X associate a A;- vertex 
graph G [20] as follows. Each vertex represents one 7^ € X, and two vertices representing 7j 
and 7j are connected by a line if they intersect (when identified with subsets of the lattice A as 
described above). An example is shown in Fig. [3l^a). A polymer occurring with multiplicity 
m, (1 <m<k\\VLX contributes m vertices pairwise connected with lines (intersects itself). 
Thus the set Pfc of all disjoint fc-polymers in (J2.33P is the set of all X whose associated graphs 
consist of k disconnected vertices (no lines). The set of all connected fc-polymers, i.e. those 
X such that no 7^ € X does not intersect at least one other 7^ G X, will be denoted by C^. 
Thus the set C^ consists of all those clusters X whose associated graphs are path-connected 
/c- vertex graphs. E.g., the cluster in Fig. [3|^a) belongs to C3. Then one has ([E], [IS]) 



In 



00 



(2.34) 



fc=i 



xeCfc 



76X 



where the index q{X) of the connected cluster X is given by 



(2.35) 



Gc oaX 



In ()2.35p the sum is over the graph and all connected subgraphs Gc on X, and 1{G(^ denotes 
the number of lines in Gc- Thus, e.g., the cluster in Fig. [3]has index q{X^ = —1 + 3x1 = 2. 



^1 ^2 ^3 



A, 



Figure 2: The set {x € 7}, in this example {xi,X2,X3}, belonging to a given connected 
diagram 7 (cf. text). 



-• • 



(a) 



(b) 



Figure 3: (a) The graph G{X) of a cluster X = {71,72,73} consisting of three mutually 
intersecting polymers, (b) The set of all connected (proper) subgraphs on X. 



Expectations (j2.4|) are correspondingly given by 



,^) = E ^ E 9(^) n C(7) . (2.36) 

k=i ' xeCklo] -yex 

In (j2.36p Cyt-[0] denotes the set of all connected /c-clusters containing (k — 1) polymers with 
activities (j2.3ip and one polymer 7© having non-empty intersection with O (the support of 
O) with activity given by 



c{ia)= j n ^^- n ^^^= ^ n ^; 



bs 



xeou7c) 63670 
(j2.36p is easily obtained from (|2.34p by writing 



ha,l,a£'^0 



o 



^^-:^lnZ^[(l + AO)]|,^,, 



where Z\[0] denotes the partition function with the insertion of the operator O. 



(2.37) 



(2.38) 



2.3 Convergence 

Having formulated our expansion as a polymer expansion for the logarithm of the partition 
function ()2.34p or for observable expectations ()2.36p one may proceed to examine its conver- 



10 



gence by a straightforward application of known results |14j , [15] concerning the convergence 
of such expansions. As it is well-known, for translation invariant systems, the following con- 
vergence criterion holds [21]: polymer expansions (j2.34p and (j2.36p converge absolutely and 
uniformly if 



Q = supV|C(7)|el^l < 1. 



79a: 



(2.39) 



In (j2.39p the sum is over all polymers 7 holding fixed a site x G 7. 

From (|2.23p - (j2.24p the propagator color- flavor matrix elements are bounded by 

2 



|C(r-r',e(x))|< 



[1 - e-M] 



-H\T — T 



(2.40) 



with |t — r'l taken as the "periodic" (shortest) distance between the two points in the timelike 
direction on the lattice torus. From (|2.28p one has 



n—uLv^L 



-fj,L 



2u 



< |DetC7-^(e(x))| <2-''V^-^ 



1 + e 



-fiL 



2u 



(2.41) 



Consider a polymer 7 build out of a set of I7I spacelike bonds with bonds of multiplicity 
m being counted m times. (Note that in fact m > 2 because of the mod Nc equal number 
of Ub and U^^s integration constraint above.) Now, Ub^, U^ are bounded by unity in color 
space. There are at most v choices for connecting a f/' or '(/' at a site (x, r) on the boundary 
of one bond to a V' or ■;/' at a boundary site (x, r') of another bond in the set via a propagator 
C{t — t'). There are I7I propagators in the diagram. (Each propagator connects two bonds 
but each bond has two boundary sites; also, note that only sites separated by an odd number 
of bonds in the timelike direction can be so connected.) From (|2.3ip - (|2.32p . (|2.13p and (|2.4Up . 
(j2.4ip then one has 



IC(7)I < (^)'l"l 



1 + e- 



-fiL 



1 



-fiL 



2u\x\ 



1 



-fiL 



xe7(T,rj) 



■^Fi 



(2.42) 



To next sum over all such polymers made of I7I spacelike bonds and attached to a fixed 
site we sum over all possible timelike separations between the bonds and over all possible 
configurations of the bonds. Now 



L/2 
r=l 



L-1 



-HT\ 



< e 



-M 



i + E' 



-/^rl 



-/^ 



r=l 



1 



-fiL 



-fJ- 



whereas the number of possible spacelike bond configurations is easily seen to be bounded by 
(2d)^l'''l. Combining with (|2.42p and noting that |x| < |7|/2 one finally has 



Q<Y Kl^l(e'^-1)-I^le 



l7l 



l7l=l 



eK 
el" -1-eK' 



where 



K = 2u'^{2df 



1 + e- 



-jiL' 



(2.43) 



(2.44) 



11 



The convergence condition Q < 1 then gives 2eK < (e^ — !)• We have thus arrived at the 
following result. For any spatial dimension d > 1 the polymer expansions \2.34^ and \2. 3(^) 
converge absolutely and uniformly in \Ks\ if 



In particular, one has: 



4ez/2(2d)2 < [tanh(^/2a^r)]'^ {e^^ - 1) . 



ln[l + 4ez^2(2d)2] <^ for r = 0. 



(2.45) 



(2.46) 



A sufficient condition on /i satifying (12.45^ for any temperature (2 < L < oo) is in fact given 
by 

In [1 + Aev'^{2df (coth(ln[l + Aev"^ {2df])y] < fi . (2.47) 




^ 



Figure 4: Striped region of convergence (chiral symmetry) given by (j2.47p in text. The 
complete chirally symmetric region is expected to include the shaded area connecting to the 
region of low /^t at high T. 



It should be noted that ()2.45p is a substantial overestimate since it was arrived at by 
ignoring the spacelike bond multiplicity constraints as well as all factors of inverse powers of 
Nc resulting from the spacelike gauge field integrations and the restrictions to odd timelike 
separations. The bound improvements from taking these into account, however, do not 
change the qualitative picture represented by Fig. HI More importantly, one may expect 
that the high fx convergence region connects to the high T convergence region [12] (through 
the shaded area in Fig. [4|). To establish this one first needs to reintroduce arbitrary {at /as) 
in the spacelike action. This ratio was conveniently set equal to one above but it provides 
crucial convergence factors at large T and small /i. After this trivial step, one needs to bound 
the gauge field dependence (at small fi) less crudely than in (I2.40p and ()2.4ip . 



3 Chiral symmetry 



An immediate consequence of the convergence of the cluster expansion (I2.36P is the existence 
of chiral symmetry within the region of its convergence. 
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Indeed, the expectation of any local chirally non-invariant fermion operator 0{x), such 
as, e.g., ijj{x)ijj{x), vanishes identically term by term in the expansion (|2.36p by the chiral 

invariance of the measure d^x- Correlation functions (0{x)0{y)) then receive vanishing 
contributions from any polymers consisting of two disjoint clusters connected to sites x and 
y, respectively. Nonvanishing contributions arise only from polymers intersecting both sites 
X and y. Now, the total number tlt of timelike bonds in propagators in any minimal such 
polymer is easily seen to be bounded from below by the minimal distance between the two 
sites. Also, there are I7I propagators in a polymer 7, and I7I < ht since each propagator is 
at least one timelike bond long. From (|2.42p then the activity of such a minimal polymer is 
bounded from above by 




-firiT 



-^ 



1 _ e-M / 2e(2(i)2 



riT 



< 



2e{2df 



riT 



where we used the convergence condition 2eK < {e^ — 1) from above. From this it follows 
that within the convergence radius of the expansion 2-point correlations are absolutely and 
uniformly bounded at any temperature from above by 



0{x)0{y) 



< Const. 



1 



\x-y\ 



2e(2d)2 
where \x — y\ is the minimum number of bonds connecting the two sites, i.e., 



lim lim 

\x—y\^fOD |As|— >c 



Oix)Oiy) 



0. 



(3.1) 



(3.2) 



In the same manner, all higher correlation functions of any chirally non-invariant fermion op- 
erators ()2.5p exhibit exponential clustering for large separations, i.e., there is no spontaneous 
breaking of the global chiral symmetry. 

4 Conclusion 

We have shown that for sufficiently large quark chemical potential strongly coupled SU{Nc) 
lattice gauge theories with Nj massless staggered fermions and at any temperature are in 
a chirally symmetric phase. This agrees with and generalizes previous results obtained by 
numerical or mean field techniques in particular cases as reviewed in section 1. 

There are several directions in which this work may be further pursued. The lower bounds 
on the magnitude of the critical chemical potential given in section 2.3 are very non-optimal 
and can be improved as outlined. As also mentioned there, it should be possible to extend 
the region of convergence in the fJ- — T plane to high T and small fx, thus connecting with the 
high T proof in [12j . This would require doing a better job of bounding timelike gauge field 
integrations at small //, and would appear quite feasible, especially in particular cases such 
as Nr = 2, 3. 
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Such improvements may be viewed within the general framework of employing the ex- 
pansion as a computation tool. In this paper we were mostly concerned with the establishing 
the convergence of the expansion rather than its computational use. Explicit computation of 
graphs involves integration over the gauge fields on spacelike bonds in the hopping expansion 
and integration over the timelike gauge field dependence in the propagators. The former, 
familiar from ordinary hopping and strong coupling expansions on the lattice, can be carried 
out systematically to quite high order. The latter will in general require numerical evalua- 
tion. Truncating the series to leading approximation will give, as already mentioned, results 
quite close to those in [9] - |llj . The expansion then provides systematic corrections to such 
leading approximations. As in the case of ordinary hopping expansions, one might hope that 
resummation of (infinite) classes of diagrams may also allow extension beyond the radius of 
convergence of the original series, and/or improved mean field or other approximations in a 
wider regime. 

A physically interesting and obvious question is whether the present treatment can be 
extended from the strong coupling limit to the finite gauge coupling region. For odd number 
of colors Nc this would appear rather problematic. For even Nc, however, and in particular for 
Nc = 2, the real positive nature of the fermionic determinants (cf. (|2.29p ) makes it possible 
to "repackage" the above expansion making use of this positivity to extend into the finite 
coupling region. We hope to explore this elsewhere. 

The author would like to thank P. de Forcrand for correspondence. 



A Appendix 

In this Appendix we derive the explicit expressions (I2.23p . ()2.24p and ()2.28p . In the following 
we set 7o = 1 in ()2.2ip for staggered fermions. We always take even L. 

Given a meromorphic function g{z) which is analytic on \z\ = 1, summation over the 
Brillouin zone p.22p can be performed by using the easily verified representation 

(both contours taken counterclockwise). Now if the integrand tends to zero as \z\~'^ or faster 
as l^l — )• oo one may deform the outer contour to infinite radius taking into account any 
poles of g{z)/z encountered in this deformation. The integral over the inner contour may be 
performed by residues of its enclosed poles. 

We use (jA.ip to evaluate the propagator (j2.20p - (j2.2ip . which, written out separately for 
T > and T < 0, has the form 

C{r, e(x)).,,, = ^^8.,6.,- J^^^ ,^^j,^Oa{.)/L^ ^^] '°' " = ^>"" ' ^^'^^ 
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For r > and fi > one then obtains 



■i 



z^ 1 



Now, since |t| < L, the integrand is sufficiently convergent for the outer contour to be 
deformed to infinite radius where it gives zero contribution. No poles are encountered in this 
deformation, the integrand having poles at z± = ±e~ "'"^Z e"'^ inside the inner contour. 
Evaluation of the inner integral by residues gives then (|2.23|) . 

For T < and // > one has from (IA.2P 

/ z''^! 1 

|2| = l+e J\z\ = l-e [Z 



.2 gi294x)/Lg2Mj [1 + ^L] 



The inner contour integral now encloses no singularities and gives zero contribution. The 
outer contour can again be deformed to infinite radius, now, however, encountering the 
integrand singularities at z± = ±e^^"-^^''^e^. Evaluating their contribution results into (j2.24p . 
The two cases are reversed for /i < (antiquark chemical potential). 

An alternative method of summation over the Brillouin zone is by use of the Fourier series 
representation of Ylm=-oo ^(^ ~ 2?TZ7r): 

L/2 



(A.5) 



1 ^ e*'='"-g(A:„,) = ^ ^ C dk 5{k - k^) e'''^ g{k) 

kmeBZ m=-(L/2)+l '^^'^ 

= ^ / dk 6{kL - {2m - 1)tt) e'^^ g{k) 

m.=-oo -^"^ 
n=-oo •^~^ 
n=-oo U-TT 

Applying (|A.5P to ([2:20]) one obtains 

CXD 

c(t, e(x)) = J2 (-1)" ^(^ + "^' ® W) (A-6) 

n=— oo 

with C(r, 0(x)) given by (I2.25|) . i..e, by the propagator (I2.20p in the L — t- oo limit of the BZ. 
For r > and /^ > 0, letting z = e*^, (j2.25p may be represented in the form: 

C(r, e(x))„,, = -5..5.,1 e--»(^)/-e- £^^^ d. ^^, _ ^_.X(x)/.^-2.j • (^-7) 

15 



Evaluating the contour integral gives: 

C'(r,e(x)),,,bi=5afe5ii[l-(-ir]e-^^'^W^/^e-^^, r > 0, m>0. (Ai 

For T = — |r| < and /i > one obtains 

.1-1 



C{T,e{^))ai,bj 






dz 



|,|=1 ^^2_g.2e.(x)/Lg2M] 



0. 



r < 0, ^i > . 



(A.9) 
(A.IO) 



This is as expected since, at fixed x, (j2.25p represents the propagator of a 1-dimensional 
static (no kinetic energy term) fermion in the background of a constant potential. For /x < 0, 
(jA.SP and (jA.lOp are reversed (antifermion propagation). 

Since L > \t\, (lAlUjl reduces ^Al\i to (ITWll . i.e., 



C(t, e(x)) = ^(-1)"C(t + nL, e(x)) . 



(A.ll) 



n=0 



(The n = term is actually absent for r < 0). Inserting the explicit form (jA.SP in (jA.lip 
then and carrying out the sum for r = ib|T| results in (j2.23p and (j2.24p . respectively. 

Integrating out the fermions in ()2.1ip . i.e., along a one-dimensional timelike chain at fixed 
x with action ()2.19p gives 



DetC-^(e(x)) 
Using the identity 



n n ^sin(^^^I^^ + eJL-^, 

a=lm=-L/2+l ^ 



L 



Nf 



L-1 

n 

n=0 



2n7r \ (-1)^/2 
sm I -^ \- z ] = „/ ^ (1 - cos(L2;)) 



L 



()A.12p becomes 

DetC~^(e(x)) = iL'^(^-l)Lu/22~Lu 



2L- 



N,, 



(A.12) 



(A.13) 



.a=l 



Nf 



vL u^L 



2-'^^e 



Nc 



.a=l 



Nf 



(A.14) 



(j2.18p was used in obtaining ()A.14p which is (j2.28p of the main text. 
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